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sj ■ ABSTRACT 



We present estimates of the photometric redshifts, stellar masses and star forma- 
tion histories of sources in the SCUBA HAlf Degree Extragalactic Survey (SHADES). 
This paper describes the 60 SCUBA sources detected in the Lockman Hole covering 
an area of ~ 320 arcmin 2 . Using photometry spanning the B band to 8^m, we find 
that the average SCUBA source forms a significant fraction of its stars in an early 
period of star formation and that most of the remainder forms in a shorter more in- 
tense burst around the redshift it is observed. This trend does not vary significantly 
with source redshift but the exact ratio of early to late mass is quite sensitive to the 
way extinction is treated in the modelling. However, the sources show a clear increase 
in stellar mass with redshift, consistent with downsizing. In terms of SED types, only 
two out of the 51 sources we have obtained photometric redshifts for are best fit by a 
quasar-like spectral energy distribution, with approximately 80 per cent of the sources 
being best fit with late-type spectra (Sc, Im and starburst). By including photometry 
at 850/im, we conclude that the average SCUBA source is forming stars at a rate 
somewhere between 6 and 30 times the rate implied from the rest-frame optical in 
a dust obscured burst and that this burst creates 15-65 per cent of the total stellar 
mass. Using a simplistic calculation, we estimate from the average star formation his- 
tory that between one in five and one in 15 bright (L» + 2mag < L op ti C ai < L* — lmag) 
galaxies in the field over the interval < z < 3 will at some point in their lifetime 
experience a similar energetic dusty burst of star formation. Finally, we compute the 
evolution of the star formation rate density and find it peaks around z ~ 2. 
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1 INTRODUCTION 

The SCUBA HAlf Degree Extragalactic Survey (SHADES) 
is a wide area extragalactic sub-millimetre (submm) survey 
conducted with S CUBA (the Submilli metre Common User 
Bolometer Array; iHolland et al. 1119991 ). The motivation for 
SHAD ES is discussed at len gth in the survey definition pa- 
per bv lMortier et a.1. 1 (|2005h . The survey comprises two sep- 
arate fields of approximately equal area, one in the Lock- 
man Hole and one in the Subaru/XMM-Newton Deep Field 
(SXDF). Up to the time of decommissioning of SCUBA in 
2005, SHADES had acquired approximately 40 per cent of 
the target area to the proposed depth of 2mJy, culminating 
in the detection of a total of 120 robu st SCUBA sources over 
- 650 arcmin 2 (|Coppin et al. 1120061 ). 

SHADES satisfies a long awaited demand for a large, ho- 
mogeneous sample of submm sources with multi- wavelength 
follow-u p data. Since their detection in the first deep submm 
surveys (ISmaii Ivison fc Blain I Il997l ; iHughes et aTTI 1 19981 ; 
iBarger et al. Ill998h . much has been learnt about the dusty 
high-redshift sources revealed by SCUBA. However, several 
major questions regarding this elusive population remain 
unanswered. 

Arguably the most important question is the relation- 
ship between SCUBA sources and present-day galaxies. Sev- 
eral clues point toward a strong link with massive low- 
redshift ellip ticals, such as their similar comoving num- 
ber d ensities (|Scott et al. 1120021 ; iDunne. Eales fc Edmunds! 
20031) . clustering properties fe.g.. lAlmaini et al. I 1200*3 : 
Blain et al. 1 120041 ) and their typically very high star for- 
mation rates which enable the rapid formation of a large 
stellar system. If this link is valid, then an immediate ques- 
tion that arises is at what stage of this transformation do we 
observe the object as a SCUBA galaxy? It is possible that 
there is more than one answer if multiple routes exist to the 
same type of final massive elliptical. For example, the aver- 
age elliptical's entire stellar population may either form in 
a single large burst or in a series of smaller bursts triggered 
by mergers. 

Crucial evidence can be provided by the star formation 
history of the average SCUBA source. However, in order 
to establish a star formation history, several ingredients are 
needed: 1) the source must be identified from its somewhat 
imprecise SCUBA position, 2) multi-wavelength data (ide- 
ally covering optical to submm) must be acquired 3) the 
source's redshift m ust be known. Deep radio surveys (e.g., 
Ilvison et ail [2002) detect somewhere between a half and 
three-quarters of SCUBA sources to give precise positional 
information. In addition, SCUBA sources can be efficiently 
identified with relatively short exposures using the Spitzer 
Space Tel escope (Spitzer) as demonstrated by s everal au- 
thors (e.g. lHuang et al. ll2004l ; lEgami et al. 1120041 ). Once ac- 
curate positions have been obtained, optical spectroscopy 
can then be carried out. 

Unfortunately, this procedure has several selection ef- 
fects. The requirement that a source be detected at radio 
wavelengths can lead to a lack of sources at high redshifts 
(z J> 3) where the radio flux falls below the detection limit 
(see for example IChapman et ail 120051 ) . There is also a 
paucity of sources at redshifts where no bright spectral fea- 
tures fall within the optical waveband, particularly over the 
redshift interval 1.2 < z < 1.9 (the 'redshift desert'). Pho- 



tometric redshifts, albeit less precise, do not suffer from the 
latter selection effect. Another advantage is that contami- 
nating flux from near-neighbours and blended sources (as 
an appreciable fraction of SCUBA sources appear to be) is 
more readily quantified in image data unlike spectroscopic 
data where small uncertainties in the slit placement can lead 
to ambiguities in deblending. 

In this paper, we make use of photometric redshifts to 
investigate the photometric properties, stellar masses and 
star formation histories of SCUBA sources in SHADES. 
Counterparts are identified through either deep radio data 
or Spitzer 24/im images. Having the option of a 24/im iden- 
tification means that the SCUBA sources considered in this 
work are not entirely subject to the strong radio selection 
function (although it may be contended that the process 
of identifying counterparts in the radio is more physically 
motivated) . 

This paper is se venth in a series of pa pers arising from 
SHADES. Paper I bv lMortier et al. I (|2005l ) de scribes the sci- 
ence goals, motivation and strategy. Paper II (|Coppin et al. I 
2006) presents the maps, catalogues an d source counts and 
describes the data reduction. Paper III (jlvison et al. II2007T I 
details the radio follow-up of the SHADES areas and identi- 
fies the radio and 2 4/im Spitzer counterpa rts to the SHADES 
sources. Paper IV l|Aretxaga et al. 1120071 ) derives photomet- 
ric redshifts of the SHADES sources using radio, submm 
and far infrared (far-IR) data. Paper V l|Takagi et al. 1120071 ) 
concerns the submm propert ies of near-IR select ed galax- 
ies in the SXDF. Paper VI l|Coppin et al. 1120071 ) presents 
350/im observations of of a subset of SHADES sources . 
Papers VII (this paper) and VIII (|Clements et al. I I2007T ) 
form a pair split by survey area: paper VII considers pho- 
tometric redshifts, stellar masses and star formation histo- 
ries of sources in the Lockman Hole whereas paper VIII 
is concerned with photometric redshifts in the SXDF. Pa- 
pers VII and VIII are divided primarily due to data propri- 
ety, but also because of different optical and near-IR cov- 
erage between the two areas. In ter ms of forthcoming pa- 
pers, paper IX (Ser jeant et al. Il2007l ) will investigate stack- 
ing of the SHADES data to statistically determine prop- 
erti es of the very faint sourc es hidden in the noise. Paper 
X l|van Kampen et al. I [2007) will study the clustering of 
SHADES sources. Finally, several further SHADES publica- 
tions are anticipated concerning 1.1mm data acquired with 
the AzTEC instrument on the James Clerk Maxwell Tele- 
scope. 

The layout of this paper is as follows. Section [2] out- 
lines the acquisition and reduction of data used in this work. 
Section discusses the results of our photometric redshift 
analysis and properties of the SEDs. In Section [4] we de- 
scribe our method of obtaining stellar masses together with 
the resulting masses, star formation rates and evolution of 
the star formation rate density. We conclude the main pa- 
per sections with a summary and discussion in Section [5] 
We provide three appendices: Appendix [X] contains multi- 
wavelength postage stamp images for each of our sources as 
well as their best fit SEDs, Appendix [B] provides descrip- 
tions for a selection of noteworthy sources and Appendix [Cl 
lists the photometry for all sources. 

Throughout this paper, we assume the following cos- 
mology: H =70kms" 1 Mpc" 1 , fi m = 0.3, Ov = 0.7. 
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2 DATA 

2.1 SHADES catalogue & counterparts 

The sources investigated in this paper were extracted from 
850/im SCUBA observations of the Lockman Hole centred 
on RA=10 h 52 m 26.7 s , Dec=57°24'12.6" (J2000). The Lock- 
man Hole data cover an area of ~ 320 arcmin 2 to an RMS 
noise level of ~ 2mJy. To ensure a robust list of SCUBA 
sources, the map reduction and catalogue generation was 
carried out by four independent groups within the SHADES 
consortium. Only sources with a signal to noise ratio of ^ 3 
(before deboosting) in at least two reductions were retained. 
This left a total of 60 sources in the Lockman Hole field with 
a probability of < 5 per cent of being spurious. For more spe- 
cific details on the reduction of the SCUBA map and sou rce 
extraction, the reader is referred to lCoppin et al. ] l|2006h . 

Follow-up 1.4GHz imaging with the Very Large Ar- 
ray and 24^m im aging with Spitzer are described in 
llvison et al. l|2007l ). Potential radio and 24/xm counterparts 
to the SCUBA sources were searched for within a radius 
of 8" and the significance of each match quantified using 
the method of lDownes et al.l (|l986l ). This method gives the 
probability, P, of a counterpart being associated with the 
SCUBA position by chance based on the separation and the 
number counts. In this way, a counterpart either in the ra- 
dio or at 24/xm is defined as being robust if its distance from 
the SCUBA position is less than 8" and P < 0.05. Non- 
robust identifications are defined as those with P > 0.05 
within 8" or those that lie within an extended search radius 
8" < r < 12.5" in the ra dio or 8" < r < 15" at 24/Ltm. 

llvison et al. I (|2007l ) lists several counterparts for many 
of the SHADES sources. In the present work, we base our 
analysis on the single, most likely counterpart for each 
source. We define these 'primary counterparts' as those with 
the lowest value of P in the radio or at 24/xm (often both). 
Although not included in any of our later analyses, for com- 
pleteness, we also compute photometric redshifts for sec- 
ondary counterparts. A secondary counterpart is defined as 
having a robust 24/xm and/or robust radio identification but 
with a numerically higher P than the primary. These are 
listed in Table [2] alongside the primary counterparts. In all 
cases, the co-ordinates listed in T able [2] are eit h er th e ra- 
dio or 24/tm co-ordinates given in llvison et al.l (|2007f ) de- 
pending on which identification has the lowest P. (There is 
one exception: Lock850.036 for which we give the SCUBA 
850^m centroid as this source has no radio or 24/im coun- 
terparts). 

2.2 Optical and near-IR photometry 

Our optical i mages in B, R, I and z were obtained with 
SuprimeCam l|Mivazaki et al. 1 120021 ) on the Subaru tele- 
scope in January 2006. The images fully cover the Lockman 
area observed by SCUBA and reach a 5a point source sensi- 
tivity of 26.8, 25.8, 25.7 and 25.0 mag (AB) in B, R, I and 
z respectively as measured in a 3" diameter aperture. To- 
tal exposure times for B, R, I and z are respectively 7200s, 
3360s, 4730s and 4800s. The seeing in the images varies from 
0.66" to 0.84" between bands with a mean of 0.76". 

Our K band image was obtained with the Wide Field 
Camera (WFCAM) on the United Kingdom Infra-red Tele- 
scope (UKIRT). The data were taken as part of the 



Deep Extra-galactic Survey, one of the five projects com- 
prising the UKIRT D eep Infra-red Sky Survey (UKIDSS 
lLawrence et al. Il2006l ). The image is a mosaic o f four sepa- 
rate quadrants observed in array number 1 (see iDve et al. I 
120061 . for further details). The exposure time of the four 
quadrants varies between 9180s and 11460s, reaching an 
average 5a point s ource sensitivity of 22.9 mag (AB) 
l| Warren et al. 1120071 ). 

We used SExtractor (jBertin fc Arnouts I [l996h to ex- 
tract sources from the optical and K band data. Only ob- 
jects with five or more interconnecting pixels lying above a 
threshold signal-to-noise of 2cr were extracted. Fluxes were 
computed within a 3" diameter aperture. 



2.3 Spitzer photometry 

Spitzefs Infra-red Array Camera (IRACi lFazio et al. 1120041 ) 
was used to obtain images at wavelengths 3.6, 4.5, 5.8 and 
8/mi. The integration time for each image was 500s, reach- 
ing 5a point source detection limits of 1.3, 2.7, 18 and 
22/tJy for 3.6, 4.5, 5.8 and 8/im respectively. The data were 
reduced using the Sp itzer Science Centre's (SSC) pipeline 
|Gordon et al. 1120051 ). 

The pipeline produces the total flux of each source. To 
match the total flux of a given source in each of the IRAC 
bands to the optical and near-IR 3" photometry accounting 
for the different PSFs, we adopted the following procedure: 
1) Fit a 2D Gaussian to the source in the IRAC image, 2) 
Scale the fitted Gaussian to the size it would have been if it 
had been observed with the SuprimeCam PSF, 3) Compute 
the flux within a 3" diameter aperture centred on the scaled 
Gaussian using the total flux output by the IRAC pipeline. 
The average of all corrections across all bands and source^] 
was +0.34 mag with a la scatter of 0.15 mag. To account for 
possible systematics introduced by this scheme, we added an 
error of 0.15 mag in quadrature to the error computed by 
the SSC pipeline for all IRAC photometry. 

In addition to IRAC imaging of the Lockman Hole, data 
were ta ken using the Mu lti-band Photometer for Spitzer 
(MIPS: iRieke et al. Iliooj ). Although we do not use MIPS 
fluxes directly in the present work, our identification of coun- 
terparts to the SCUBA sources reli es on 24/xm MIPS d etec- 
tions (see previous section and also llvison et al."ll2007h . 



2.4 Cross matching procedure 

To obtain the multi-wavelength list for all the SHADES 
counterparts, the sources on the different images were po- 
sition coincidence matched using a tolerance of 1". All 
matches were carefully verified by eye to correct for obvious 
mismatches, spurious detections and blended photometry. 
In cases where blended photometry was identified, we ad- 
justed the the deblending threshold and cleaning parameter 



1 Since the correction applies to the total flux which is practically 
insensitive to the PSF and since sources are effectively brought 
to the same seeing in each IRAC waveband, the scatter in the 
correction between wavebands for a given source is primarily due 
to the image noise and wavelength dependent source morphology, 
not the varying PSF between wavebands. 
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in SExtractor on a source by source basis until deblended 
photometry was obtained. 



3 PHOTOMETRIC REDSHIFTS AND 
SPECTRA 

3.1 Method of determination 

We obtained photometric redshifts for the SHADES sources 
by a pplying the HyperZ redshift code l|Bolzonella et alTl 
2000) to our nine band photometry (B, R, I, z, K, 3.6/im, 
4.5/^m, 5.8/xm and 8/zm) described in the previous section. 
The template galaxy spectral energy distributions (SEDs) 
package d with HyperZ are bas ed on the local SEDs mea- 
sured bv lColeman et al~l (|l980h . with an extrapolation into 
the infrared using the results of spectral synthesis models. 
The lack of any empirical basis for these templates at wave- 
lengths > lfim is clearly unsatisfactory since our photom- 
etry includes both near-infrared and mid-infrared measure- 
ments. Therefore , we c onstructed our own set of templates. 
iMannucci et al. I (|200lh list empirical SEDs over the wave- 
length range 0.1/im to 2.4/im for the Hubble types E, SO, 
Sa, Sb and Sc. We extended these out to 10/xm using the 
average SED s for disk g alaxie s and elliptical galaxies listed 
in Table 3 of iLu et al. I l|2003h . 

One disadvantage of these templates is the lack of a tem- 
plate for an irregular galaxy, and we therefore retained the 
HyperZ Im template, extending this into the mid- i nfrare d 
using the average disk-galaxy SED from ILu et al~l l|2003h . 
In terms of the variation in SED shape, particularly the av- 
erage slope between 0.2/im< A < 2fim, the jump from the Sc 
template to the Im template is larger than the progression 
through the earlier type SEDs. Therefore, we introduced a 
seventh intermediate template with a composition of 60 per 
cent Sc and 40 per cent Im to bridge this jump. 

Finally, we introduced a starburst (SB) and quasar 
(QSO) t emplate. For th e SB te mplate, we extended the spec- 
trum of iKinnev et al. 1 (1 19961) to lo n ger a nd shorter wave- 
lengths using a iBruzual fc Chariot I (|2003l ) spectrum corre- 
sponding to a O.lGyr duration starburst as observed at the 
end of that time. Since this extension lacks the dusty spec- 
tral featur es seen in the m id-IR, we added these from the 
spect ra of ILu et al. I (|2003 | ). T he QSO template was taken 
from iBrotherton et al. I 1 200 if ) a nd extended longward of 
7500 Ausing the QSO spectrum of lElvis et al~l i| 19941) . This 
is an average SED and therefore does not cater for a variable 
power-law slope. Similarly, we chose not to include templates 
containing a mix of QSO and starburst activity. Our pref- 
erence is for a more clear-cut approach: The QSO template 
that we use is sufficiently different from the other templates 
to be well suited for picking out obvious QSO candidates. 

The nine SED templates are plotted in Figure Q] 
Throughout this paper, we refer to SED templates by their 
number; 1=E, 2=S0, 3=Sa, 4=Sb, 5=Sc, 6=Sc+lm, 7=Im, 
8=Starburst, 9=QSO. 

To determine the redshifts, we allowed HyperZ to search 
over the range < z < 6. We used a Calzetti reddening law 
to account for dust extinction, allowing the rest-frame V- 
band attenuation, Ay, to vary from to 5 mag in steps of 
0.1 mag (requiring re-compilation of HyperZ with a larger 
array). No cut-off was imposed on the absolute K band mag- 
nitude (AB) computed by HyperZ (see Section [4| although 
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Figure 1. Template SEDs used with HyperZ for determination 
of the photometric redshifts. The template Sc+Im is a hybrid 
SED composed of 60 per cent Sc and 40 per cent Im to bridge the 
jump between the Sc and Im types (see text). SED templates are 
referred to throughout this paper by their number; 1=E, 2=S0, 
3=Sa, 4=Sb, 5=Sc, 6=Sc+Im, 7=Im, 8=Starburst, 9=QSO. 



we imposed a minimum magnitude error of 0.05 mag for all 
filters to account for zero point errors. We configured Hy- 
perZ to treat non-detected sources as having zero flux with 
a la error equal to the flux sensitivity of the corresponding 
filter. (We also tried a more stringent configuration whereby 
the error is set to half the sensitivity with negligible effect 
on the results) . All filter transmission profiles were modified 
where appropriate to account for the wavelength-dependent 
quantum efficiency of the detector, the mirror reflectivity 
and atmospheric attenuation. 

We note that we include the quantity p(z) in Table [2] 
and its corresponding p(z2) for secondary redshifts. This 
quantity is the integral of the area under the peak of the 
redshift probability distribution normalised by the total area 
within < z < 6 for a given source. The redshift probabil- 
ity distribution is computed from the \ 2 curve output by 
HyperZ. The most likely redshift (i.e., lowest \ 2 ) i s taken 
as the primary solution, not the peak with the largest en- 
veloped area. Figure [5] shows p(z) versus z for two extreme 
cases in Table [2] At one extreme is Lock850.003 where a 
very robust photometric redshift with a small error and no 
secondary solution is obtained and at the other extreme is 
Lock850.100 where three almost equally likely solutions ex- 
ist. We find that the majority of sources have a single well 
defined peak. 



3.2 Photometric redshifts and SED results 

3.2.1 Photometric redshifts 

Table [2] lists the photometric redshifts obtained for 51 pri- 
mary counterparts (and 12 secondary counterparts) of the 
60 SHADES sources and the lower plot in Figure [3] shows 
their distribution. The best fitting SEDs are shown in Figure 
IA1I The median redshift of the sample is z — 1.52 in close 
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Figure 2. Two extreme examples of the photometric redshift 
probability distribution of the SHADES sources studied in this 
paper. The top panel corresponds to the source Lock850.003 
which has a very robust photometric redshift. Conversely, the 
bottom panel corresponds to Lock850.100 which has three almost 
equally likely photometric redshifts. The grey shaded area indi- 
cates the 90% confidence range in primary z as given in Table 

m 



agreement with the median p hotometric redsh i ft me asured 
in the SHADES SXDF field bv lClements et al. I i|2007h . Red- 
shifts for nine of the sources could not be determined due 
to insufficient photometric points (either because they are 
too faint and therefore not detected at certain wavelengths - 
we stipulate a minimum of three fluxes per source to deter- 
mine a redshift - or because they are heavily blended with 
a neighbouring source). 

A simple measure of the appropriateness of our SED 
templates to the photometry is the distribution of x 2 - This 
analysis depends on the accuracy of our photometric errors 
although for simple aperture fluxes, these are reliably deter- 
mined. The distribution for all our SED fits is perfectly con- 
sistent with a x 2 distribution of 5 degrees of freedom. This 
is as it should be for nine photometric constraints (including 
upper limits) and four fitted parameters (z, Av, SED nor- 
malisation and SED type) implying that our SED template 
set provides a suitable match to the SCUBA sources. 

We have compared the photometric redshifts output by 
HyperZ with three alternative determinations (see Table [T] 
for a summary). Firstly, still applying HyperZ bu t using 
the synthetic spectra of iBruzual fc Chariot I l|2003t ) pack- 
aged with it, we found generally good agreement although 
a slightly larger spread i n the dis t ribut ion of x 2 - Secondly, 
we used the BPZ code of lBenitezl | |200Cf ) with our own SED 
templates. Although this code does not allow for attenu- 
ation by dust, it uses a Bayesian approach incorporating a 
luminosity prior. We set the prior to the luminosity function 
obtained from the Hubble Deep Field North and found ex- 
cellent agreement with the results of HyperZ when no dust 
attenuation was included. However, the inclusion of dust at- 
tenuation in HyperZ gives significantly different redshifts for 
a large fraction of the sources and values of x 2 that are a 
factor of two times lower on average. As expected, this shows 
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Figure 3. Bottom panel: Histogram of photometric redshifts com- 
puted for the 51 primary SHADES sources listed in Table[2] Mid- 
dle panel: Variation of SED type with redshift. Top panel: Varia- 
tion of attenuation by dust with redshift. 



that dust has a significant effect on the SEDs of the SCUBA 
sources (also reflected in the values of Av seen in Table [2}. 

In the third and final case, we compared the HyperZ 
redshift estimates with our own photometric redshifts ob- 
tained using the star formation history analysis presented 
in Section 14.11 Figure H] shows the comparison. The agree- 
ment is very good with an average |Az|/(l + z) of 0.09 (for 
comparison, the average a z /(l + z) for the HyperZ redshifts 
is 0.06). Since the two methods are completely independent, 
with one using empirical SED template fitting and the other 
synthetic spectra constructed from a best fit star formation 
history, this supports the reliability of our redshifts. 

Unfortunately, there are very few SCUBA sources in the 
Lockman Hole with robust spectroscopic redshifts against 
which we can compare our photometric redshifts. Never- 
theless, in Figure [5] we show a comparison with the 7 
that currently exist. The only irreconcilable discrepancy is 
Lock850.017 with z phot = 1.06± and z spcc = 2.24. The 
optical counterpart to this source coincides exactly with the 
radio and spectroscopic position, is not blended with any 
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Photo-z type 


< |Az| > 


< cr z > 


HyperZ 




0.10 


HyperZ+BC 


0.34 


0.14 


BPZ 


0.21 


0.31 


SFH 


0.18 


0.16 


submm 


0.78 


0.60 



Table 1. Comparison of the different photometric rcdshifts 
(photo-zs) considered in Section 13.2.11 < <r z > is the median 
1(7 uncertainty in the photo-zs of all 60 SHADES sources and 
< |Az| > is the median absolute difference between the photo- 
zs listed in Table [2] and the others discussed. The photo-z types 
listed are: HyperZ - the photo-zs in Table [2] i.e., obtained us- 
ing HyperZ with the empirical SED templates discussed in Sec- 
tion 13.11 HyperZ + BC - again using HyperZ but with the 
iBruzual fc Chariot I j2003ft SED templates; BPZ - using the code 
bv lBenitez I j2000h with the emperical SED templates of Section 
13.11 SFH — our own photo-zs from the star f ormation history 
analys is of Section l4.ll submm - the photo-zs of Arctxaga ct al. 
d2007h . 



Figure 4. Comparison of the photometric redshifts computed 
by HyperZ and those from the star formation history analysis of 
Section l4.ll (ltr error bars plotted). 



Lock850.014 
Lock850.03.' 



* 2 



Lock850.004b' 




Lock850J)l 
Lock850.016 



-Lock850.017 



spec 



Figure 5. Comparison of photometr i c reds hifts and robust spec- 
trosco pic redshifts from llvison et al. ] (20051) and lChapman et al. I 
(2005). Error bars show the 90 per cent confidence range on the 
photometric redshifts. 



neighbours and is detected with a very high signal to noise 
at all wavelengths. The photometry is extremely well fit by 
the Im+Sc template with a reduced x 2 of 0.28 ± 0.10 and a 
single sharp peak in the probability distribution of the fit. 
However, the optical morphology shows a relatively diffuse 
tail extending ~ 3" to the south-west of the nucleus of this 
source. A possible cause of the discrepancy could therefore 
be that there are two closely aligned objects at different 
redshifts. Our photometry is dominated by the compact nu- 
cleus which could be a foreground source, whereas the spec- 



troscopic redshift could be based on emission lines from the 
more diffuse background source. 

We draw particular attention to the interesting source 
Lock850.041b. The radio co-ordinates of this source corre- 
spond to the brightest but slightly less favoura ble of two ra- 
dio so urces within 8" of the SCUBA centroid. Ilvison et al. I 
(2005) were unable to measure a spectroscopic redshift at 
the radio co-ordinates of Lock850.041b but did measure a 
spectroscopic redshift of z = 0.69 for an elliptical galaxy 
lying ~ 2" to the south-west of the radio co-ordinates. We 
determine a photometric redshift of z — 2.22 for the ra- 
dio source and find that the photometry is well fit by a 
QSO SED template. Combined with the optical morphol- 
ogy which exhibits arc-like structure (the arcs being a result 
of a relatively bright extended host galaxy) , these facts sug- 
gest that Lock850.041b is very likely a stron g ly len sed QSO. 
This supports the claim of I Chapman et ai"~l l|2002l ) that the 
source is a plausible lens. 

Figure [6] compares our redshift s with the photomet- 
ric redshifts of lAretxaga et al.~l (|2007l ) determined using the 
850/im/1.4GHz spectral index (their 'z^ot')- We find a large 
but not significant offset between both sets of redshifts (not 
including sources with lower limits) . Comparing the sources 
with z su bmm derived from more than two photometric points 
(the filled circles in Figure [6}, our redshifts are lower by 
Az — 0.37 ± 0.39 on average, where the error here is the 
error on the mean. Extending this comparison to include all 
sources in Figure [6] we find that our redshifts are lower by 
an average of Az = 0.42 ± 0.35 (or alternatively, as given in 
Tabled] the median value of \Az\ is 0.78). 

This offset manifests itself more strongly when compar- 
ing the pe ak of our total r edshift distribution at z ~ 1.5 
to that of Aretxaga et al. I j2007) which lies at z ~ 2.3 
(Chapman et al.l 20051 , also find that the peak of their 
spectroscopic redshift distribution lies at z ~ 2.3). The ori- 
gin of this discrepancy is not entirely clear, although the 
offset is reduced if the comparison is limited to the red- 
shift distributions of only those sources plotted in Figure 
(i.e., those with redshifts determined by both the opti- 
cal and submm/radio methods). This suggests that part 
of the discrepancy is a selection effect. If this is the case, 
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submm 

Figure 6. Comparison of photometric redshifts computed in this 
work (,z n hnt) with the photometric redshifts of lAretxaga et al. I 
1120071 ) derived using the 850/xm and 1.4GHz spectral index 
(^ubmm)- Open squares correspond to redshifts derived solely 
from 850^tm and 1.4GHz photometry whereas filled circles corre- 
spond to redshifts derived with additional far-IR/submm photo- 
metric data. Both vertical and horizontal error bars denote the 
90 per cent confidence range. The tabbed ends of the arrows show 
the lower limits set by lAretxaga et al. I [120071 ). 



then it is interesting to n ote how the two di s tribut ions 
compare to the models of Ivan Kampen et al~l (|2005l ). A 
KS test shows that our redshift distribution is consistent 
with the massive merger model whereby all SCUBA galax- 
ies are the result of a violent merger of two galaxy sized ha- 
los. C onversely, the redshift distribution of lAretxaga et ahl 
(|2007h is best fit by the phenomenological model, which in- 
cludes a mixtur e of starbursts due to m ergers and quiescent 
star formation. lAretxaga et al. I l|2007h note that their red- 
shifts are even better matched by the distribut i on pr edicted 
from the semi-analytical model of lSilva et al. I l|2005t ) which 
jointly describes the formation and evolution of spheroids 
and QSOs. An implication of this might therefore be that 
SCUBA sources more readily detected in the optical are 
more likely to be the result of massive mergers, whereas a 
pure submm selected sample is more likely to include sources 
with quiescent star formation. 



3.2.2 SED properties 

Figure IA1I in Appendix [X] shows the best fit SEDs to the 
SHADES sources for which we have been able to estimate 
redshifts. The SED type and Av given in Tableware plotted 
in Figure [3] as a function of redshift. There are two distinct 
trends. Firstly, there is a clear evolution in SED type, with 
late types being found predominantly at higher redshifts and 
early types at lower redshifts. Secondly, sources at z > 2.5 
have significantly lower attenuation than sources at z < 2.5: 
Dividing the top panel of Figure [3] by the line Av = 1, it can 
be seen that sources at z < 2.5 are split equally, whereas at 
z > 2.5, 11 sources have an attenuation Av < 1 with none 



15 



10 



123456789 
SED type 

Figure 7. Histogram of best fit SED types for the 51 primary 
SHADES sources listed in Table [2] SED templates are; 1=E, 
2=S0, 3=Sa, 4=Sb, 5=Sc, 6=Sc+Im, 7=Im, 8=SB, 9=QSO. 



having Av > 1. We have verified that this is not a result 
of the way HyperZ deals with shorter waveband dropouts 
by reproducing this trend independently of the configured 
method (see Section |3,1|I , 

Such a decline in attenuation at high redshifts is consis- 
tent with the rate of dusty star formation having peaked at 
an earlier epoch (see Section [43]), although there is a strong 
selection effect due to the sensitivity of the data. This de- 
cline is also reflected in the declining number of SHADES 
sources found beyond z > 2.5. This is not entirely a result 
of the usual selection effects since our redshifts are not sub- 
ject to the optical spectroscopic redshift desert discussed in 
Section [T] nor are they totally reliant on there being a radio 
detection since the ID procedure also relies on detections 
at 24/^m. A similar decl ine in attenuation due to dust has 
recently been found by iBuat et"aL Their findings 

indicate a monotonic decrease in dust attenuation in lumi- 
nous infrared galaxies from z = out to z = 2 for a fixed 
star formation rate (see also lXu et al. Il2007h . 

Figure [7] shows the distribution of SED type (listed in 
Table[2]) for the 51 primary sources. As the figure shows, the 
SHADES sources are dominated by late type SEDs. Over 
80 per cent of sources have SEDs of type Sc or later and 
approximately two thirds of these are best fit with either 
the Im or SB template. There are only two sources that are 
best fit with a QSO template: Lock850.075 and Lock850.081 
(there is also the secondary Lock850.041b). This is entirely 
consistent with the low fraction of AGN dominated S CUBA 

20021: 



Ivison et al. 



sources found in previous studies (e.g., 
Waskett et al~ll2003l : lAlmaini et aLl2003l : lAlexander et"al"1 
2005t l 



4 STELLAR MASS 

4.1 Method of determination 

To determine stellar ma sses of the SHADES source s, we use 
the synthetic spectra of iBruzual fc Chariot I |200of) . Rather 
than assume an a priori model for the star formation his- 
tory (SFH) of each galaxy, we divide the history of each 
galaxy into blocks of fixed fractional time (see Section |4,2[) 
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Table 2. Photometric redshifts of the Lockman Hole SHADES sources. Those sources with an ID suffix 'b' are the secondary counterparts 
listed here for completeness but discounted from all analysis (see Section l2.1|l . Codes in column labelled 'C indicate whether the source 
has a radio and/or 24^tm counterpart. First character corresponds to radio and second to 2Afim with 'R' indicating a robust identification, 
'N' a non-robust id entification and '-' no detection (see Section |2.1> . All co-ordinates are either the radio or 24/^m co-ordinates from 
llvison et al. (2007) depending on which is the more likely identification (exception is Lock850.036 - see text). Columns z + and z give 
the 90 per cent uncertainty range on the photometric redshift, z. Column 'SED' lists best fit template number (see text). Absolute AB 
magnitude listed under Mjf. x 2 of the SED fit and the fractional area contained under the peak of the probability distribution for x 2 > 
p(z), are given in columns If and 12 (see text). Colum ns 13 to 15 lis t secon dary r edshift solutions and th eir corresponding x 2 an d P- 
Last column gives spectroscopic redshifts from either (I) llvison et al. I ||2005|) or (2) Chapman et al. ] ||2005| ) (robust redshifts in bold). 



ID 


C 


RA 




Dec 








z 


z 


z+ 


SED 


A v 


M K 


x 2 W 


P(*) 


22 


X 2 (*2) 


p(z2) 


^spec 


Lock850.00f 


RR 


10 


52 


01.249 


+57 


24 


45.76 


4 


.21 


2.41 


4.45 


6 


0.2 


-24.15 


0.76 





.87 


- 




- 


2.15 1 


Lock850.002 


RR 


10 


52 


57.084 


+57 


21 


02.82 


3 


.23 


2.76 


3.34 


7 


0.7 


-24.10 


0.22 





.98 


- 


- 


- 


- 


Lock850.003 


RR 


10 


52 


38.299 


+57 


24 


35.76 


1 


.21 


1.12 


1.25 


1 


0.0 


-22.91 


0.69 


1. 


.00 


- 


- 


- 


3.04 1 


Lock850.003b 


RR 


10 


52 


38.401 


+57 


24 


39.50 


1 


.51 


1.24 


1.61 


2 


0.1 


-23.36 


0.99 


0. 


.36 


0.71 


1.02 


0.37 


- 


Lock850.004 


RR 


10 


52 


04.079 


+57 


26 


58.52 


2 


.66 


2.25 


2.85 


7 


0.9 


-24.71 


0.26 





.73 


3.16 


0.73 


0.26 


0.53 1 


Lock850.004b 


RN 


10 


52 


04.226 


+57 


26 


55.46 


2 


.11 


1.59 


2.50 


6 


0.9 


-24.26 


0.67 


1. 


.00 


■ 


- 


- 


1.48 1 


Lock850.005 


-N 


10 


53 


02.696 


+57 


18 


21.95 


1 


.20 


0.87 


1.24 


2 


0.1 


-22.96 


0.52 


0. 


.96 


■ 


- 


- 


- 


Lock850.006 


RR 


10 


52 


04.013 


+57 


25 


24.20 


2 


.01 


1.81 


3.18 


8 


2.0 


-23.67 


0.60 





.98 


- 






- 


Lock850.007 


RR 


10 


53 


00.956 


+57 


25 


52.06 




_ 


_ 


- 


- 


_ 


_ 


_ 














Lock850.008 


-K 


10 


51 


53.690 


+57 


18 


34.90 




- 


- 


■ 


■ 


- 


- 


_ 






■ 


- 


- 




Lock850.009 


RR 


10 


52 


15.636 


+57 


25 


04.26 


1 


.62 


1.54 


1.87 


8 


2.1 


-23.88 


0.75 


0. 


.97 


■ 


- 


- 


1.85 1 


Lock850.009b 


-R 


10 


52 


15.730 


+57 


25 


01.70 





.91 


0.82 


1.01 


6 


0.4 


-23.01 


0.43 


1 


.00 










Lock850.010 


R- 


10 


52 


48.992 


+57 


32 


56.26 




- 




- 


- 


- 


- 






. 


- 








Lock850.010b 


-N 


10 


52 


48.270 


+57 


32 


51.00 


1 


.33 


■ 


1.39 


7 


1.1 


-22.09 


■ 





.98 


■ 


- 


- 


- 


Lock850.0II 


-N 


10 


51 


29.160 


+57 


24 


06.80 


1 


.33 


1.23 


1.42 


7 


1.6 


-21.99 


0.32 


0. 


.61 


1.75 


0.58 


0.36 


■ 


Lock850.012 


RR 


10 


52 


27.579 


+57 


25 


12.46 


2 


.03 


1.83 


2.33 


7 


1.5 


-24.07 


0.75 





.93 


2.77 


1.56 


0.07 


2.14 1 


Lock850.013 


-N 


10 


51 


31.770 


+57 


31 


41.20 





.48 


0.31 


0.55 


7 


0.1 


-21.70 


0.61 


1 


.00 


_ 


- 


- 


- 


Lock850.014 


NN 


10 


52 


30.717 


+57 


22 


09.56 


2 


.76 


2.47 


3.08 


6 


0.1 


-25.17 


0.45 


0. 


.81 


1.97 


0.89 


0.19 


2.61 1 


Lock850.015 


RR 


10 


53 


19.271 


+57 


21 


08.45 


2 


.53 


2.08 


2.94 


5 


0.0 


-24.24 


1.45 





.99 


- 


- 


- 


- 


Lock850.015b 


R- 


10 


53 


19.025 


+57 


21 


09.47 




_ 


_ 


_ 


_ 


_ 


_ 


_ 




_ 










Lock850.016 


RR 


10 


51 


51.690 


+57 


26 


36.09 


1 


.25 


1.05 


1.38 


5 


1.2 


-24.26 


0.65 





.88 


0.71 


1.31 


0.12 


1.15 1 


Lock850.017 


RR 


10 


51 


58.018 


+57 


18 


00.27 


1 


.12 


1.00 


1.38 


6 


1.3 


-23.39 


0.28 


1. 


.00 


- 


- 


- 


2.24 1 


Lock850.0I7b 


-R 


10 


51 


58.480 


+57 


18 


01.20 





.36 


0.30 


0.44 


8 


2.1 


-17.91 


3.55 


0. 


.69 


0.15 


3.98 


0.29 


■ 


Lock850.018 


R- 


10 


52 


27.778 


+57 


22 


18.18 




_ 


_ 


_ 




_ 


_ 


_ 




_ 


_ 


_ 


_ 


1.96 1 


Lock850.019 


-R 


10 


52 


36.090 


+57 


31 


19.60 


3 


.69 


3.56 


3.75 


' 


0.3 


-24.35 


3.56 





.88 


3.00 


4.30 


0.11 




Lock850.02I 


-R 


10 


52 


56.790 


+57 


30 


37.90 





.91 


0.43 


3.24 


1 


1.9 


-21.13 


0.08 


0. 


.90 


0.22 


0.77 


0.07 




Lock850.022 


-R 


10 


51 


37.090 


+57 


33 


16.90 


2 


.68 


2.55 


2.90 


8 


0.1 


-24.11 


0.89 


0. 


.93 


1.86 


2.24 


0.07 


- 


Lock850.023 


-N 


10 


52 


14.976 


+57 


31 


53.62 





.57 


0.53 


0.68 


1 


0.3 


-23.28 


1.67 


1 


.00 


_ 


_ 


- 




Lock850.024 


RR 


10 


52 


00.445 


+57 


20 


40.16 


1 


.40 


1.25 


1.49 


8 


1.5 


-24.12 


0.72 


1 


.00 










Lock850.026 


RR 


10 


52 


40.698 


+57 


23 


09.96 


2 


.86 


1.84 


3.19 


7 


1.0 


-23.43 


1.11 





.98 


- 


- 


- 




Lock850.027 


-N 


10 


52 


03.450 


+57 


18 


19.30 


1 


.12 


1.03 


1.20 


7 


0.4 


-22.47 


1.17 


1 


.00 


■ 


- 


- 


- 


Lock850.028 


NN 


10 


52 


57.667 


+57 


30 


58.71 


1 


.11 


1.08 


1.20 


8 


2.1 


-23.29 


0.35 


1 


.00 






- 




Lock850.029 


N- 


10 


51 


31.305 


+57 


20 


40.28 




_ 


_ 


- 


_ 


_ 


- 


- 








- 


- 




Lock850.030 


RR 


10 


52 


07.490 


+57 


19 


04.01 


2 


.06 


1.85 


2.56 


6 


0.2 


-23.01 


0.74 


0. 


.82 


0.52 


1.47 


0.17 


2.69 2 


Lock850.03I 


RR 


10 


52 


15.989 


+57 


16 


19.34 


2 


.12 


1.95 


2.31 


6 


1.0 


-24.44 


0.81 


1. 


.00 










Lock850.033 


RN 


10 


51 


55.470 


+57 


23 


12.77 


2 


.58 


1.98 


2.98 


6 


0.3 


-23.35 


0.38 





.55 


1.08 


0.83 


0.33 


2.66 1 


Lock850.034 


RN 


10 


52 


14.202 


+57 


33 


28.30 


3 


.42 


3.15 


3.69 


7 


0.6 


-24.05 


0.92 


1 


.00 










Lock850.035 


N- 


10 


52 


46.655 


+57 


20 


52.54 






























Lock850.035b 


-N 


10 


52 


45.940 


+57 


20 


51.40 


1 


.59 


1.41 


1.71 


7 


0.6 


-22.68 


1.03 


1. 


.00 










Lock850.036 




10 


52 


09.335 


+57 


18 


06.78 






























Lock850.037 


R- 


10 


51 


24.342 


+57 


23 


36.18 






























Lock850.037b 


NR 


10 


51 


24.595 


+57 


23 


31.08 


1 


.51 


1.46 


1.66 


8 


1.7 


-23.61 


1.01 


1 


.00 










Lock850.038 


RR 


10 


53 


07.060 


+57 


24 


31.60 


1 


.29 


1.17 


1.39 


8 


1.2 


-23.67 


0.85 


1. 


.00 










Lock850.039 




10 


52 


25.505 


+57 


16 


08.54 






























Lock850.040 


RN 


10 


52 


01.721 


+57 


19 


17.00 


2 


.02 


1.87 


2.14 


6 


0.6 


-23.88 


2.92 





.98 










Lock850.041 


RR 


10 


51 


59.760 


+57 


24 


24.94 


1 


.01 


0.91 


1.11 


8 


1.0 


-24.13 


2.02 


1 


.00 










Lock850.041b 


RR 


10 


52 


00.248 


+57 


24 


21.69 


2 


.22 


2.05 


2.44 


9 


0.9 


-24.37 


0.86 


0. 


.52 


0.32 


1.17 


0.26 




Lock850.043 


NR 


10 


52 


56.561 


+57 


23 


52.80 


1 


.72 


1.67 


1.91 


7 


1.4 


-23.49 


0.99 


1 


.00 










Lock850.043b 


NR 


10 


52 


56.576 


+57 


23 


58.62 


1 


.14 


1.04 


1.42 


8 


1 + 


-23.96 


0.39 


1 


.00 










Lock850.047 


-N 


10 


52 


34.850 


+57 


25 


04.60 





.37 


0.33 


0.42 


1 


0.0 


-20.21 


4.49 





.91 


3.97 


4.80 


0.09 




Lock850.048 


-R 


10 


52 


56.030 


+57 


32 


42.30 





.31 


0.24 


0.37 


8 


1.6 


-21.28 


0.81 


1. 


.00 










Lock850.052 


RR 


10 


52 


45.808 


+57 


31 


19.86 


1 


.21 


1.15 


1.30 


1 


0.4 


-24.00 


0.66 


0. 


.57 


0.96 


0.96 


0.43 




Lock850.052b 


-R 


10 


52 


46.160 


+57 


31 


20.20 


1 


.40 


1.28 


1.60 


5 


1.5 


-23.62 


0.81 





.84 


0.71 


1.63 


0.12 




Lock850.053 


-R 


10 


52 


40.290 


+57 


19 


24.40 


1 


.55 


1.45 


1.74 


7 


0.5 


-23.19 


1.22 


1 


.00 










Lock850.060 


-N 


10 


51 


43.900 


+57 


24 


43.60 


1 


.92 


1.12 


2.34 


1 


0.1 


-22.38 


1.11 


0. 


.60 


2.69 


1.37 


0.16 
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Table 2 — continued 



ID 




C 


RA 




Dec 






z 


z 


4- 


SED 




M K 


2 t \ 
X («) 


P(«) 


z2 


2 / o\ 

X (22) 


p(z2) Z spcc 


Lock850 


063 


RR 


10 


51 


54.261 


+57 


25 


02.55 


4.73 


4.62 


4.89 


7 


0.2 


-25.50 


0.14 


1 


00 




- 






Lock850 


064 


NN 


10 


52 


52.320 


+57 


32 


33.00 


1.19 


1.11 


1.36 


8 


1.1 


-22.45 


0.29 


1 


00 




- 






Lock850 


066 


-N 


10 


51 


39.570 


+57 


20 


27.10 


1.41 


1.01 


1.94 


8 


0.3 


-20.70 


0.05 





54 


2 


15 


0.27 


0.33 


Lock850 


067 


-B 


10 


52 


08.870 


+57 


23 


56.30 


1.46 


1.33 


1.80 


1 


0.0 


-21.89 


0.71 


1 


00 




- 






Lock850 


070 


NN 


i n 
1U 


01 


47.894 


+57 


30 


44.37 


0.53 


0.46 


0.62 


7 


1.2 


-20.86 


0.84 


1 


00 










Lock850 


071 


RN 


10 


52 


19.086 


+57 


18 


57.87 


1.91 


1.76 


2.07 


6 


0.3 


-23.18 


1.44 


1 


00 






- 


- 


Lock850 


073 


RR 


10 


51 


41.992 


+57 


22 
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and determine the amount of star formation in each block. 
We choose five blocks as a compromise; too few gives an un- 
acceptably low SFH resolution whereas too many introduces 
high degeneracy between blocks and reduces the number of 
sources for which a stellar mass can be computed due to 
missing photometry. Blocks are approximately logarithmi- 
cally sized to account for the fact that a galaxy's SED is 
more strongly influenced by more recent star formation ac- 
tivity. 

We use the 199 4 Padova stellar evolutionary tracks 
Bertelli et al. I [1994) with a Salpeter initial mass function 



55) 



Salpeter 1 1 19551 ). Starting with a simple stellar population 



(SSP) SED, Lf SF , we generate a composite stellar popu- 
lation (CSP) SED, LI, for the ith block of constant star 
formation in a given galaxy using 



1 



dt'Lf p (r(z)-t') 



(1) 



where the block lasts from time ti-i to ti in the galaxy's 
history and r is the age of the galaxy (i.e., the age of the 
Universe today minus the lookback time to the galaxy). The 
SSP SED is normalised to one solar mass hence the constant 
Ati = ti — ti-i ensures that the CSP is also normalised to 
one solar mass. We then redden the CSP SED of each block 
using the extinction Ay via the relationship 



-0Ak(X)A v /R v 



(2) 



where fc(A) is given by the Calzetti law for starbursts 
ijCalzetti et al. 1120001 ), 



fc(A) = 



2.659(-2.156 + 
2.659(-1.857 + 



1.5U9 _ 0.198 , 0.011 \ j 

(for 0.12/xm < A < 0.63^m) 

1 -j ± ) + Rv 

(for 0.63itm < A < 2.2/xm) 



(3) 



with Rv = 4.05 and A in microns. To match our template 
SED coverage, we assume that the longer wavelength half 
of the function applies up to 10/^m, and we linearly extrap- 



olate the shorter wavelength half down to 0.01/im using the 
average slope between 0.12/^m and 0. 

The flux (i.e., photon count) that would be observed in 
filter j from a given block i at the redshift z of the galaxy is 



Fl] ~ 47Td? 



dA 



ALUA/(l + z))T J (A) 
(1 + z) he 



(4) 



where is the luminosity distance and Tj is the transmis- 
sion curve in filter j (this includes telescope and atmospheric 
throughput as well as detector response). We then find the 
normalisation a; of each block of star formation by minimis- 
ing the x 2 function 



3 



(5) 



where F° hB is the flux observed in filter j from the galaxy and 
(Tj is its error. We treat non-detections in the same way as we 
configured HyperZ to, i.e., the flux is set to zero and assigned 
a la error equal to the sensitivity of the corresponding filter. 

We use a downhill simplex method to minimise \ 2 ■ 
To prevent non-physical solutions we impose the constraint 
ai > 0;i = l,iVbiock- Since the CSP SED from which Fij is 
computed is normalised to one solar mass, the quantity ai 
is the amount of stellar mass in solar units formed by the 
galaxy in the time interval Ati. The total stellar mass of the 
galaxy is therefore simply 



JVblock 



(6) 



The error on M* is calculated by summing in quadrature 
the error on each ai derived from the fit. 

In minimising \ 2 m equation we have the choice 
of either minimising only the quantities ai and holding Ay 

2 We have investigated the effect of different extrapolations on 
our results and find negligible dependence. Since the cut-off wave- 
length of 0.12^tm does not reach the central wavelength of our 
shortest filter (B band) until z ~ 3, this is not surprising. 
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Figure 8. Variation of stellar mass of the SHADES sources 
with rest- frame absolute K band magnitude (AB). Mass error 
bars are determined by a Monte Carlo analysis and include un- 
certainty due to metallicity (see text) and errors on Mr de- 
rived from redshift and photometric uncertainty. The outlying 
point corresponds to Lock850.081, best fit with the QSO tem- 
plate. Excluding this source, the straight line fit is log 10 M* = 
(0.59 ± 0.05) - (0.468 ± 0.004) M K . 



Figure 9. Variation of stellar mass of the SHADES sources with 
redshift. Mass error bars are determined by a Monte Carlo anal- 
ysis and include uncertainty due to metallicity (see text). The 
continuous curve shows the mass limit transformed from the K 
band sensitivity limit using equation J7) and our starburst SED 
template. The lack of objects within 7xl0 n < M*/M Q < 4xl0 12 
at z < 2 compared to the 12 objects within 2 < z < 5 is consistent 
with downsizing. 



and 2 fixed at the values determined by HyperZ, or minimis- 
ing all parameters a%, Av and z. Although the latter option 
seems the more compelling given its self-consistency, as we 
show in Figure [4] the redshifts obtained in both cases show 
very good agreement on the whole (see Section [3. 2. Nev- 
ertheless, there are some significant differences and since the 
templates used with HyperZ are empirical and include dust 
features, we choose the former option for all analysis in this 
paper. 

4.2 Stellar mass results 

We determined the stellar mass of each SHADES source 
using the prescription given in the preceding section. An as- 
sumption of this is that the metallicity is fixed as the source 
evolves. Also, the same metallicity is used for all sources. 
To account for these limitations, we treated the metallicity 
as a source of uncertainty, repeating the calculation of stel- 
lar mass for two extremes . To s et this range, we used the 
results of ISwinbank et al~l (|2004h who measured the metal- 
licity of 30 high redshift SCUBA sources. They found an 
approximately solar metallicity on average with a la scat- 
ter of ~ 0.25dex. We therefore performed two calculations 
of stellar mass for each source, one using the spectra of 
iBruzual fc Chariot I l|2003l ) with metallicity Z = 0.4 Z and 
a second with Z = 2.5 Zq. The resulting difference in mass 
was taken as the 68 per cent confidence range and combined 
with the Monte Carlo error as described below. 

Figure [5] shows the stellar masses computed for the 
SHADES sources, plotted against the rest-frame absolute 
K band magnitude. The correlation is very strong with lit- 
tle scatter. Of the two sources best fit with QSO templates, 
Lock850.081 is an outlier (as labelled) but Lock850.075 lies 
within the main trend with a slightly higher than median 
stellar mass of 3.9 x 10 11 M®. The most plausible explana- 
tion for this difference is that unlike Lock850.081, the pho- 
tometry of Lock850.075 is not dominated by central QSO 



emission. If we omit Lock850.081, linear regression gives the 
following relationship 

log 10 M„ = (0.59 ± 0.05) - (0.468 ± 0.004)M K . (7) 

The rest-frame K band flux is dominated by old stars and 
therefore, as has been appreciated for some time, provides a 
good mass estimator for such populations. The small scatter 
seen in Figure [8] shows that rest-frame K also provides a 
good calibrator for the total stellar mass in SCUBA galaxies. 
This is a direct consequence of the fact that the average 
SCUBA galaxy appears to have an early and a late block of 
star formation which makes approximately equal amounts 
of old and new stars (see next section). 

The average stellar mass of our sample is 1O 11 ' 8±O1 M0. 
This compares favourably with an average of 10 11A±0 ' i Mp, 
from the study of 13 SCUBA galaxie s by iBorvs et al. I 
(2005). Similarly. [Swinbank et alT(|2006l ) estimated dynam- 
ical masses of eight submm galaxies, finding an average of 
1O U - 7±O - 3 M . In terms of mass-to-light (M/L), the straight 
line fit expressed by equation (0 gives M/L cx L 017±0 01 , 
indicating that the mass increases more rapidly than the 
luminosity. 

In Figure [9] we show how stellar mass varies with pho- 
tometric redshift. To obtain the mass error, we performed 
a Monte Carlo analysis with 1000 realisations using solar 
metallicity SEDs. In each realisation, we computed the mass, 
randomly sampling the fluxes and redshift for each source 
using the measured errors and assuming a normal distribu- 
tion. The lcr scatter in stellar mass for each source was then 
added in quadrature to the lcr error resulting from the un- 
known metallicity as described above. These two errors are 
typically approximately equal. 

The continuous curve in Figure [9] shows the mass sensi- 
tivity limit. We estimate this by firstly computing the abso- 
lute rest-frame K band magnitude corresponding to the 5<r 
point source sensitivity of 22.9 mag (AB) at each redshift as- 
suming our Im template SED. This is then transformed into 
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stellar mass using the fitted relationship given in equation 
0. 

The horizontal dashed lines in Figure [9] show the arbi- 
trary selection 7 x 10 11 < M*/M Q < 4 x 10 12 . These limits 
select the 12 most massive sources. All 12 sources lie within 
2 < 2 < 5. Since the comoving volume over this redshift 
interval is nearly twice that over < z < 2 (and since our 
mass sensitivity barely affects the 2 < z < 5 selection), one 
would expect from simple Poisson statistics to find ~ 6 ± 2 
sources within < z < 2. However, there are no SHADES 
sources with masses greater than 7 x 10 1 Mq at z < 2. This 
is evidence in favour of downsizing, whereby star formation 
in the Universe progressively shift s to smaller systems at 
later times (e.g.. lCowie et al.lll996h . 



4.3 Star Formation Rate 

In deriving these stellar masses, we divided the history of 
each galaxy into the periods (0-0.45)r, (0.45-0. 70)r, (0.70- 
0.85)r, (0.85-0.95)r and (0.95-l)r where r is the age of the 
galaxy, taken as the age of the Universe today minus the 
lookback time to the source. All sources were therefore as- 
sumed to start forming immediately after the Big Bang. 

The grey-shaded histogram in Figure [TO] shows the nor- 
malised star formation ratqj (SFR) in the five blocks, aver- 
aged over all 51 primary SHADES sources with photometric 
redshifts. As in the determination of stellar mass, the er- 
ror, shown in the figure by the upper and lower dashed his- 
tograms, incorporates the uncertainty due to the unknown 
metallicity and the scatter from the Monte Carlo analysis. 
Clearly, the SFR is on average dominated by a short burst 
close to the epoch at which the source is observed. Since this 
result is derived from mainly rest-frame optical photometry, 
the SFR in the last block will be suppressed due to obscura- 
tion by dust. We estimate the effect of this in Section [47372] 

Determination of the average SFR in this way leads to a 
surprising result. Even though the SFR is dominated by the 
later stages in the average source's history, the quantity of 
stellar mass created, i.e., the area under each block in Figure 
[TOl is split approximately evenly between these later stages 
and the earliest stage (considering only the optical+IR de- 
rived results for now). Specifically, during the earliest pe- 
riod, (42 ± 12) per cent of the total stellar mass was formed, 
compared to (39 ± 6) per cent during the fourth period. In 
the context of our model (see below), this therefore indi- 
cates that the average SHADES source has al ready formed 
a sign ificant fraction of its stars (also noted bv lBorvs et al. I 
2005) and that it is undergoing a second major episode of 
star formation at the epoch at which it is observed. 

To verify the robustness of this result, we have carried 
out two tests. In the first test, we simply repeated the cal- 
culation twice but using two different sub-sets of SHADES 
sources. In the first sub-set, we selected those sources which 
meet the criteria x 2 { z ) ^ 1-5 an( A P( z ) ^ 0.9. This leaves 27 
out of the 51 primary sources with redshifts. In the second 
sub-set, we selected sources which have a robust identifi- 
cation in either the radio or at 24/im or both. This gives a 



3 Wc define the normalised SFR as simply the fractional stellar 
mass formed in each block divided by the block width. Section 
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14. 51 considers the absolute SFR. 
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Figure 10. The normalised SFR averaged over all 51 primary 
SHADES sources (shaded histogram, left ordinate). Upper and 
lower dashed histogram gives the 68 per cent confidence range 
allowing for uncertainty in metallicity and including the Monte 
Carlo scatter described in the text. The data points (right ordi- 
nate) give the fraction of stellar mass created in each block on 
average. The last block also shows our estimate of the SFR and 
stellar mass created when the star formation obscured at opti- 
cal/IR wavelengths is taken into account. The cross-hatched his- 
togram and open histogram are computed from the submm flux 
assuming a cold and hot SED respectively (see Section I4.3.2H . 
The bold rhombi give the corresponding stellar mass fractions. 
In this last block, both the mass fraction and SFR derived from 
the submm flux have been added to those derived from the opti- 
cal/IR. 



sub-set of 35 sources. The repeat calculation gives an almost 
identical result for both sub-sets, with ~ 40 per cent of their 
stellar mass being created during the first period and ~ 40 
per cent being created in the fourth period. 

In the second, more sophisticated test, we investigated 
how well our method reconstructs an input SFH. We gen- 
erated three synthetic source catalogues each matching the 
number of SHADES sources and each made with a differ- 
ent average SFH. We used the Bruzual and Chariot SED 
library (with Salpeter IMF and Z = Zg) to generate syn- 
thetic photometry in our nine wavebands. For each source, 
we took the redshift and Ay determined by HyperZ and as- 
signed photometric errors to the fluxes using a flux scaling 
relationship derived from the real SHADES catalogue. We 
then applied the analysis of Section 14.11 to assess how well 
the input average SFHs could be recovered. 

Figure [TT] shows the results. The three panels from top 
to bottom respectively show a constant SFH, a SFH under- 
going exponential decay and a SFH mimicking that exhib- 
ited by the SHADES sources. The error bars in this plot 
include the Monte Carlo scatter due to flux and redshift un- 
certainty as before. The errors are smaller than in Figure 
[10] since each source is assigned exactly the same SFH, un- 
like the real source population which will inevitably have an 
intrinsic scatter. In all three cases, the reconstructed SFH 
very faithfully reproduces the input SFH, indicating that 
the trend seen in Figure [10] is not an artifact of our method. 
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Figure 11. Input (dashed line) and reconstructed (continuous 
line) average SFHs. Top: Constant SFH. Middle: Exponential de- 
cay SFH. Bottom: SFH matching the average exhibited by the 
SHADES sources. The error bars include the Monte Carlo scat- 
ter due to flux and redshift uncertainty. 



4.3.1 SFR Evolution 

We have investigated the evolution of the trend seen in the 
SFR by dividing the SHADES sources equally into a low 
and high redshift bin. We defined the low redshift bin by 
1 < z 1.9 and the high redshift bin by 1.9 < z ^ 5 
(we limited this analysis to z > 1 since the eight sources 
at z < 1 poorly sample this large fraction of the Universe's 
history). Figure Q2] shows the SFHs for the two redshift bins. 
We also show for the final SFH block (i.e., 0.95-1 of the 
fractional age) our estimate of the star formation that is 
completely obscured by dust (Section l4.3.2p . In this case, for 
each redshift slice, we take the average of the SFR computed 
from the hot and cold submm SED. 

Figure [12] shows that there is little difference between 
the fractional formation rate of stellar mass at high and low 
redshift. In both the high and low redshift sample, the initial 
and late peaks in the rate of stellar mass formation persist. 
The late peak for the low redshift sample is marginally less 
pronounced than for the high redshift sample but still con- 
sistent within the errors. Of course, in absolute terms, the 
late peak in the low redshift bin corresponds to a very dif- 
ferent star formation rate to that in the high redshift bin for 
two reasons. Firstly, as Figure[9]shows, the median mass of a 
SHADES source in the high redshift bin is ~ 3 times higher 
than the median mass in the low redshift bin. Secondly, the 
proper time interval over which the late peak spans is ap- 
proximately twice that for a source at z = 1.5 than for a 
source at z = 3.5. Therefore, on average, the late peak in 
the high redshift bin corresponds to ~ 5 — 10 times the ab- 
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Figure 12. Average fraction of stellar mass formed as a fraction 
of SHADES source age for two different redshift bins (~ 22 ob- 
jects per bin). Data points for the lower redshift bin are offset 
horizontally by —0.025 for clarity. Errors account for uncertainty 
in metallicity and include the Monte Carlo scatter described in 
the text. The triangular and square points correspond respec- 
tively to the low and high redshift slices and include an estimate 
of the stars formed that are hidden by dust (see Section 14.3.21 - 
these points are averaged over the hot and cold SEDs). 

solute star formation rate of the peak in the low redshift 
bin. 

4.3.2 SFR: Dust correction 

The results presented so far are based mainly on rest-frame 
optical photometry. Although we have made a correction for 
extinction, in any system undergoing massive star forma- 
tion, some of the stars are completely obscured at optical 
wavelengths. This is a particular problem for the SHADES 
sources because we observe them as a result of a massive 
starburst. In this section, we use the observed submm fluxes 
to estimate the amount of hidden star formation. 

We assume that the far-IR/submm luminosity of a 
SHADES source is the result of a starburst extending in 
time over the period of the last SFH block and that this is 
completely hidden in the optical/near-IR wavelength range. 
We first calculate a pair of bolometric luminosities for each 
source, and L l cold , using the submm flux and a hot and 
cold SED. These SEDs are the hottest and coole s t IRA S 
galaxy SEDs from the sample of iDunne fc Ealesl (|200ll r1. 
We then calculate the hidden star-formation rate by taking 
the ratio of this luminosity and the bolometric luminosity 
-^optical f° r each source i of the CSP for the final block gen- 
erated during the analysis of the previous section: 

SFR = £[hot,cold] / -^optical • (8) 

1 i 

This gives an upper or lower limit depending on whether 
the hot or cold SED is used. Since the error caused by the 

4 We choose to estimate the range of submm bolometric flux 
in this way for homogeneity, rather than use additional observa- 
tional constraints (such as radio fluxes). The coolest SED is that 
of NGC958 dominated by cold dust at 20K and the hottest is 
IR1525+36 with a mix of dust at 26K and 57K in the ratio 15:1. 
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uncertainty in metallicity is insignificant compared to the 
range of SFRs spanned by the hot and coid SEDs, we ignore 
its contribution in this case. 

The last block in Figure [10] shows the additional hid- 
den SFR implied by the submm flux. The cross-hatched 
histogram bar is the hidden SFR assuming the cold SED 
and the empty histogram bar corresponds to the hot SED. 
Over the last 5 per cent of its history, the hidden SFR of 
the average SHADES source is somewhere between 6 and 
30 times the SFR implied by its optical flux. Similarly, the 
bold rhombi show the range in the extra amount of stellar 
mass created in this hidden burst which is 15-65 per cent 
of the total stellar mass (Note, in the figure, both the stel- 
lar mass fraction and SFR estimated from the submm have 
been added to the normalised quantities derived from the 
optical+IR.) If we therefore include the hidden star forma- 
tion, somewhere between 50 per cent and 65 per cent of the 
stellar mass was created in the last 15 per cent of the average 
SHADES source's history. 

4.4 The relationship of SCUBA sources to the 
field 

Figure [10] shows that the average SHADES source is ob- 
served during the most intense period of star formation 
the source has ever experienced. An interesting question is 
whether all galaxies endure such a phase, or whether the 
SHADES sources are rare in this respect. Based on the most 
recent two blocks in the average SHADES source's SFH, we 
can make the very crude statement that if all galaxies have 
star formation histories like the average SHADES galaxy 
found in this study, we would expect to see somewhere be- 
tween 5 per cent and 15 per cent of all galaxies undergoing 
a phase of highly energetic dusty star formation at any one 
epoch. 

To test whether this is the case, we constructed a plot 
of the apparent / band magnitude versus redshift for all the 
galaxies detected in our / band image of the Lockman Hole. 
Redshifts of the field galaxies were obtained in exactly the 
same way as the SHADES sources using HyperZ with the 
same template SEDs. We formed a master I band catalogue 
using SExtractor then coincidence matched sources at all 
other wavelengths using a radial tolerance of 1". We rejected 
stars using the CLASS-STAR parameter output by SExtrac- 
tor, only retaining objects with CLASSJ3TAR< 0.95 (leav- 
ing 87 per cent of objects). Objects in the vicinity of highly 
saturated stars were also excluded to avoid deblending prob- 
lems and pixel bleeding (predominantly in the optical) as 
were objects detected in less than five wavebands for consis- 
tency with the SHADES sample. This leaves a total of 17,000 
objects detected with 5 or more pixels above a threshold of 
2cr in the full 320 sq. arcmin SHADES area. 

Figure [l3l shows the I — z plane for the SHADES sources 
(large black points) and the full field galaxies (small grey 
points). For comparison, the continuous curve shows the 
observed I band flux of an L» galaxy computed assuming 
no evolution and using our I band filter response and Im 
template. The norm alisation of this curve is taken from 
iBlanton et ai~l |2003t ) who measure M» = 5 log h - 20.82 
in I. The plot illustrates two important facts: 1) the bright 
envelope of full field galaxies follows the L* galaxy track 
fairly closely up to and beyond the median redshift, further 
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Figure 13. Apparent I band magnitude versus redshift for the 
51 primary SHADES source counterparts listed in Table [2] (black 
points) compared to all galaxies in the field (grey points). The 
continuous curve corresponds to a non-evolving L* Im galaxy. 
The dashed curves are the Im track offset by +2 mag and -1 mag 
and the dashed straight line at Iab = 27.0 is the 90 per cent 
completeness limit. 



Redshift Fraction 

0^2<1 (1.24 ±0.51) per cent 

l<z<2 (0.72 ± 0.17) per cent 

2 < z < 3 (0.86 ± 0.26) per cent 



Table 3. Evolution of the number of SHADES sources as a frac- 
tion of galaxies in the full Lockman Hole field. Both the SHADES 
sources and full field galaxies are selected by L* + 2mag < L < 
L» — Imag and Iab 27.0. Errors assume Poisson noise only. 

demonstrating that our photometric-redshifts are at least 
reasonable, 2) the average SHADES source is significantly 
brighter than the average field galaxy. 

Before the question regarding the rarity of SCUBA 
sources can be addressed, we must make a few further con- 
siderations. The SHADES galaxies have bright / band mag- 
nitudes and are therefore among the most luminous galaxies 
at every redshift. We wish to restrict our comparison to field 
galaxies that are similarly luminous. We therefore limit both 
the SHADES and full field samples to those objects that lie 
within -1 mag and +2 mag of the L, track (see Figure [13)) . 
Secondly, we apply an upper limit of 7ab = 27 where the 90 
per cent completion limit takes effect. Finally, the majority 
of the SHADES sources are at z < 3 hence we limit our com- 
putation of the fraction to three redshift bins, ^ z < 1, 
1 < 2 < 2 and 2 sC z < 3. 

Applying these constraints, we find in all three bins that 
approximately 1 per cent of galaxies in the field is a SHADES 
source. However, since we established from the average SFH 
that only 5-15 per cent of SHADES sources would be seen 
during their most active period, we can make the very ap- 
proximate estimate that somewhere between one in five and 
one in 15 bright galaxies in the field will experience a highly 
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energetic phase of dusty star formation. This fraction shows 
little evolution over ^ z < 3 and is consistent within the 
(large) errors with being constant at all redshifts in this 
range. 

4.5 SFR density evolution 

We determined the evolution of the SFR density using the 
51 SHADES sources for which we were able to obtain pho- 
tometric redshifts. The sources were divided into six equal 
redshift bins and the co-moving volume density of the total 
SFR computed in each bin. For the SFR of each galaxy, we 
took the average SFR in the last two SFH blocks with and 
without the submm flux contribution. 

The thin dashed line in Figure [14] shows our result 
from the optical+IR photometry. The 68 per cent confi- 
dence range is indicated by the grey shading and accounts 
for the uncertainty in metallicity, the Monte Carlo scatter 
as in previous sections and Poisson noise. The thin continu- 
ous lines, in addition to the optical+IR photometry, include 
the submm flux assuming the hot and cold submm SEDs 
described in Section 14.3.21 The peak in both cases lies at 
z ~ 2, although the optical+IR peak is sharper than the 
optical+IR+submm peak, implying that the submm flux is 
emitted over a broader redshift interval. For comparison, we 
plot the SFR density m easu red from only submm flux by 
IChapman et al. I (|2005l ) and lAretxaga et al. I |2007l '). Both 
of these show a peak around 2 < z < 3, at a slightly higher 
redshift than our peak but not significantly so. The height 
of our optical+IR+submm peak is also consistent with both 
studies. 

Compared to op tical/UV observations such as the Keck 
Deep Field survey of lSawiciki fc Thompson" (120061) an d the 
Lyman break ga l axy s urveys of Steidel et al. I 1 19991 ) and 
iGiavalisco et al. I (|2004h . the SFR density is ~ 5 - 10 times 
lower for the SHADES sources. This suggests that at all 
redshifts, most of the star formation is occurring in more 
modest systems than the SHADES galaxies (which are at the 
bright end of the luminosity function). However, since the 
surface number density of SHADES sources in our sample 
is ~ 100 times lower out to a redshift of z ~ 4 than the 
galaxies observed in these other surveys, the average rate 
of star formation per SHADES source is ~ 10 — 20 times 
higher. 



5 SUMMARY AND DISCUSSION 

Using nine-band photometry ranging from B to 8[im, we 
have determined best fit SEDs, photometric redshifts and 
stellar masses as well as the average star formation history 
for 51 SCUBA sources in the Lockman Hole. We find a me- 
dian redshift of z = 1.52 (with all objects falling in the 
interval 0.22 < z < 4.73), c onsistent with the med ian pho- 
tometric redshift found by Idements et al. I (|2007T ) for the 
SHADES SXDF field. 

Approximately 80 per cent of sources are best fit with 
late type spectra ranging from Sc to starburst. Only two 
out of the 51 objects are best fit with a QSO spectrum. 
Four object s are fit with an e arly t ype (E or SO) SED, con- 
sistent with Idements et al. I (|2007T ) who find two ellipticals 
in 33 sources. This is unexpected, especially since these four 



sources all have low extinction (Ay $S 0.3). Misidentifica- 
tions aside, one possibility could be that these objects have 
a very mixed stellar population. Young star forming regions 
with a high column density of dust would be detected by 
SCUBA whilst being heavily obscured in the rest-frame op- 
tical. Meanwhile, old stellar regions with relatively little dust 
would be detected in the rest-frame optical giving the ap- 
pearance of an early type SED. 

A surprising find is that the average SCUBA source has 
already built up a significant fraction of its stellar mass in 
an early period of star formation with the majority of the 
remainder being created in a much later and more intense 
burst. This is consistent with the findings of iBorvs et al. I 
(2005) who concluded that the average SCUBA source al- 
ready has a massive population in place by z — 2.2. In- 
cluding 850^im photometry indicates that a further 15-65 
per cent of the total stellar mass is created in an ongoing 
burst of dust obscured star formation. The most recent 5- 
15 per cent of the average SCUBA source's history (220 to 
660 Myr on average) shows the highest rates of star for- 
mation ever experienced by the source. Coupled with the 
fact that ~ 1 per cent of bright field galaxies selected by 
L» + 2mag < L op ticai < L* — lmag over the redshift range 
< z < 3 appear to be SCUBA galaxies, we estimate that 
between one in five and one in 15 of these galaxies will at 
some point in their lifetime experience a similar energetic 
burst of dust obscured star formation. 

We find that the trend of an early and late formation of 
stellar mass with little intermediate activity does not differ 
between high (1.9 < z < 5) and low (1 < z < 1.9) redshifts. 
This suggests that the typical SCUBA source is a snapshot 
of a system at the same point in its history undergoing the 
same transformation process. This transformation is from a 
system with an already established, mature stellar popula- 
tion to a system with at least as much stellar mass again. 
This is an intrigui n g resu lt and fits neatly with the conclu- 
sion of iBell et af~l (|2004 ) that the average elliptical galaxy 
has doubled its mass since a redshift of z = 1. In addition 
to the early and late periods of stellar mass formation, it is 
possible that there could also have been so called 'dry merg- 
ers' between two systems containing old stars but little gas. 
In such an event, very few new stars would be formed and 
this therefore would not manifest itself in the star formation 
history at the point of the merger, but would enhance the 
stellar mass inferred from the early source history. 

There is a distinct lack of SHADES sources in the 
redshift interval < z < 2 with stellar masses greater 
than M* = 7 x 10 11 Mq , compared to 12 sources within 
2 < z < 5 above the same mass limit. This is clear evi- 
dence in favour of a downsizing scenario, where star forma- 
tion shifts to progress i vely s maller systems as the Universe 
ages. Idements et al. I (|2007T ) have found exactly this trend 
in the SXDF. 

Finally, we have determined the evolution of the star 
formation rate density using optical, IR and submm photom- 
etry. The peak occurs in the vicinity of z ~ 2, consistent with 
that d etermined from submm only studies l|Chapman et al.l 
120051 ; lAretxaga et al. I 120071 ') and that derived from op- 
tical/UV photometry (e.g., ISawiciki fc Thompson 1 120061 ) . 
Since our sample amounts to a total of only 51 sources, we 
are limited by Poisson noise throughout most of the work 
presented here. This is especially true at z < 1 where we 
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Figure 14. The SFR density of the SHADES sources. The thin dashed line is the mean density estimated from optical/near-IR 
photometry, with the grey shading giving the 68 per cent confidence range. The thin continuous lines show the same results but 
including the SFR implied from the submm flux (upper for hot SED, lower for cool SEP; see S cction l4.3.2l l. Plotted for compariso n are 
the dust-corrected S FR densities derived from optical/UV observations of Sa wiciki &: Thompson I J2006 ). Giavalisco et al. ] (|2004h and 
ISteidel et al. I lll999f). Also plotted are the SFRs determined from SCUBA galaxies by ICh*apm^net*liL* j2005Tl (extrapolated down to 
Irmly - ) and I Aretxaga et al. I <j2007T > (corrected by completing the luminosity function at 60/im). 



detect only eight sources. Future investigations with sub- 
stantially more sensitive instruments such as SCUBA-2 and 
Herschel will vastly improve this shortfall. 



APPENDIX A: SED PLOTS AND POSTAGE 
STAMPS 

Figure [All plots the best fit template SEDs to the SHADES 
sources with redshifts listed in Tabled 

Postage stamp images for all wavelengths are illustrated 
in Figure lA"2l 

APPENDIX B: NOTES ON SOURCES 

In this appendix, we provide descriptions for a selection of 
noteworthy sources. 

Lock850.001: This is a distinct source detected in RIz 
and all IRAC bands apart from 5.8/im. There is an 
unambiguous detection at 24/im and 1.4GHz. The source is 
probably blended with a fainter neighbour to the W in the 
3.6 and 4.5/im data. The SED and redshift are derived 
from deblended photometry but using blended photometry 
makes little difference to the resulting SED and redshift of 
z = 4.2. 

Lock850.003: There are two highly likely counterparts to 
this source, both with robust radio and 24/im detections. 



The primary is the fainter but more probable in both the 
radio and at 24/im with an offset of ~ 1" from the SCUBA 
position compared to ~ 3" of the secondary. The two 
counterparts have consistent photometric redshifts of 
z — 1.2li 'o9 and z = 1.51^g'2?. These are inconsistent 
with the spect roscopic redshift of z = 3.04 deemed 
non-robust bv llvison et al. I ((2005) . 

Lock850.004: There are two robust counterparts, both 
with significant detections at all wavelengths, in particular 
at 24/im and 1.4GHz. The primary counterpart has the 
lowest P in both the radio and at 24/im. The spectroscopic 
redshift for the primary is considered non-robust. A third 
source lying ~ 9" to the NNW is brighter than the primary 
in the radio and at 24/im and has a photometric redshift of 
z — 3.0 with an early type SED. 

Lock850.005: The most likely counterpart is a faint 
non-robust 24/im source lying ~ 5" from the SCUBA 
position in a SSE direction. There are no radio 
counterparts detected for this source. 
Lock850.007: There is a single robust radio and 24/im 
counterpart but all wavebands show strong blending with a 
bright neighbour which cannot be reliably deblended. This 
source is therefore omitted from the photometric redshift 
analysis. 

Lock850.008: Optical photometry of this source is 
affected by pixel bleeding from a nearby bright star and at 
IRAC wavelengths, the source is blended with a neighbour 



16 S. Dye et al. 







1- 

1 $ 


o 






\ ,_ 




Lock850.009b ~ 


d 





0.2 0.5 1 2 5 10 0.2 0.5 1 2 5 10 0.2 0.5 1 2 5 10 0.2 0.5 1 2 5 10 

A (/Mm) A (um) A (jtim) A (jl.m) 





0.2 0.5 1 2 5 10 0.2 0.5 1 2 5 10 0.2 0.5 1 2 5 10 0.2 0.5 1 2 5 10 

A (jLtm) A (/im) A (fJ,m) A (/itn) 





Lock850.017 





0.2 0.5 1 2 5 10 0.2 0.5 1 2 5 10 0.2 0.5 1 2 5 10 0.2 0.5 1 2 5 10 

A (urn) A (tim) A (iim) A (urn) 



Figure Al. Best fit SEDs for the SHADES sources determined by HyperZ. Photometric points plot 1<t error bars. Table [2] lists 
photometric redshifts, SED type, Ay and absolute K band magnitude for each source. 
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Figure A2. FIGURES OMITTED DUE TO astro-ph UPLOAD LIMIT Postage stamp images (24" X 24") for all 60 SHADES 
sources. Columns from left to right correspond to; B, R, I, z, K, 3.6/tm, 4.5/tm, 5.8/tm and 8/tm. The cross-hair shows the location of 
the primary counterpart and the diamond, where present, shows the secondary. The 10" radius circle is centred on the SCUBA position. 



of comparable brightness. The source is therefore omitted 
from the photometric redshift analysis. 
Lock850.009: The primary counterpart is the only robust 
radio source and the most likely robust 24/mi source. The 
secondary counterpart included in Table [2] has a robust 
24/im detection (P = 0.043) offset 3.8" from the SCUBA 
position but is not detected in the radio. 
Lock850.010: The robust radio counterpart reported in 
llvison et al. I ll2007h is only detected in the first two IRAC 
bands, hence no photometric redshift can be determined 
for this source. The secondary counterpart in Table [2] is the 
most likely 24/im counterpart, is detected at all 
wavelengths (apart from the radio), has an offset of ~ 8" 
SSW from the SCUBA position and is not robust with 
P=0.31. 

Lock850.011: There are no robust radio or 24/im 
counterparts for this source. The most likely counterpart is 
not detected in the radio but is detected at 24/im with a 
near-robust P=0.053. 

Lock850.015: This source has two robust radio 
counterparts, both close to each other (~ 2") and close to 
the SCUBA position (both offset by ~ 2"). The primary 
counterpart is brighter in the radio than the secondary and 
is the only source of the two to be detected at 24/im. The 
secondary is included in Table [5] by virtue of having a 
robust radio detection but is not detected at any other 
wavelength and hence has an undetermined photometric 
redshift. 

Lock850.017: At optical wavelengths, the primary 
counterpart to this source has a morphology comprising a 
centrally bright compact nucleus and a more diffuse tail 
extending to the SW. An explanation for the discrepancy 
between the photometric and spectroscopic redshifts may 
be that the tail is a separate background object; the 
photometry is dominated by flux from the compact object 
whereas emission lines seen in the spectra may originate 
from the diffuse object in the background. The secondary 
counterpart with a robust 24/im detection listed in Table [2] 
is offset ~ 2" due E, is very faint at all wavelengths and is 
best fit with the SB template but with an unacceptable x 2 - 
Lock850.021: No radio counterparts are detected for this 
source. The only robust 24/im detection is taken as the 



primary counterpart but this is only detected in the four 
fRAC wavebands. The only other 24/im source detected 
within 15" of the SCUBA position is non-robust, has an 
offset of ~ 6" SW, is brighter than the primary at all 
wavelengths and has a photometric redshift of z — 1.4 with 
an Im type SEIX 

Lock850.023: llvison et al~l (|2007T l report two non-robust 
24/im counterparts with approximately equal values of P. 
There are no radio sources detected. The slightly more 
probable (numerically lower P) 24/im source is taken as 
the primary counterpart and has a very faint 24/im flux 
and an SED that corresponds toaz~ 0.1 elliptical 
(consistent with its optical morphology). The less probable 
24/im source has a stronger 24/im flux and photometry 
consistent with a starburst at z ~ 2.7. 
Lock850.027: The primary counterpart has no radio 
detection but a non-robust faint 24/im detection 
(106 ± f5/iJy) offset ~ 6" to the N from the SCUBA 
position, consistent with an Im galaxy at z = 1.1. Another 
possible counterpart, not listed in Table [2] since it is offset 
~ 11" to the SE, has a faint non-robust 1.4GHz detection 
and a brighter but still non-robust 24/im detection 
(196 ± 13/iJy). The optical morphology and photometry of 
this source are consistent with an SO galaxy at z ~ 0.6. 
Lock850.029: No counterpart is detect ed at the pos ition 
of the most likely radio source noted by llvison et al. I 
l|2007l ) at any other wavelength. Another source lies ~ 8" 
NW from the SCUBA position and is detected at 24/im 
with P=0.15. This source is detected with high S/N at all 
wavelengths, particularly at IRAC wavelengths (but is not 
detected in the radio) and has photometry consistent with 
a Sa galaxy at z = 1.4. 

Lock850.030: This counterpart has strong radio and 
24/im flux and lies ~ 3" to the SW of the SCUBA position. 
It is faint but unambiguous at almost all other 
wavelengths, giving rise to a relatively uncertain 
photometric redshift, marginally consistent with the 
measured spectroscopic redshift. 

Lock850.035: Th e most probable cou nterpart to this 
source reported by llvison et al. I (|2007h is a non-robust 
radio source with a probability P = 0.065 offset to the SW 
by ~ 5" from the SCUBA position. This source is not 
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detected at any other wavelength, therefore we include a 
secondary counterpart in Table [2] as an alternative 
possibility, this object having the most probable 24/im 
detection seen at other wavelengths. This secondary 
counterpart is offset by ~ 9" SW of the SCUBA position, 
has irregular optical morphology and is well fit by an Im 
SED at z = 1.6. 

Lock850.036: No radio or 24/im counterpart detected for 
this source. Directly E at an offset of ~ 8" from the 
SCUBA position lies an irregular object which shows 
bright emission at IRAC wavelengths. Surprisingly, the 
photometry of this source is well fit by an elliptical SED 
with a redshift of z = 0.9. 

Lock850.037: The most probable counterpart reported by 
llvison et al. I ((2007) is detected only in the radio. Tabled 
includes a secondary counterpart with a robust 24/im 
detection and a non-robust radio detection lying ~ 6" SE 
from the SCUBA position. The secondary is the only 
source detected within 10" at 5.8 and 8/im and the 
brightest of two at 24/im. 

Lock850.041: There are two very likely counterparts, 
both with robust detections at 24/im and in the radio. 
Both have a high S/N at all wavelengths. The most likely 
counterpart has photometry consistent with az~l 
starburst, whereas the slightly less favourable counterpart 
(Lock850.041b) is very well fit by a QSO template. As 
explained in the text (see Section \3. 2. 1[1 . Lock850.041b is 
most likely a z ~ 2 .2 QSO strongly magnified by a z = 0.7 
jlvison et al. 1120051 ') elliptical lens 2" to the SW. 
Lock850.043: There are two robust 24/im counterparts 
detected at all wavelengths for this source. The primary is 
the nearer of the two with a lower P value and offset W by 
~ 5" from the SCUBA position. The counterpart has a 
significant but not quite robust radio detection and a 
weaker 24/im flux. The secondary counterpart is slightly 
blended in B with two neighbours lying to the E, becoming 
less blended at longer wavelengths where it becomes 
dominant. The secondary counterpart lies ~ 8" NW of the 
SCUBA position, has a stronger 24/im flux but a very 
similar photometric redshift of z = 1.20. 
Lock850.047: The counterpart for this source has a very 
disturbed optical morphology with a long (~ 4") 'tail' 
extending NNE. The radio and 24/im position are both 
located at the same point in the 'tail'. The photometry was 
measured in an aperture centred on the radio position and 
is probably contaminated by the brighter more compact 
source from which the tail appears to emanate. This most 
likely causes the poor SED fit and hence uncertain 
photometric redshift of z = 0.4. 

Lock850.052: Both primary and secondary counterparts 
have consistent photometric redshifts (z ~ 1.2) and both 
are very well fit with similar SED types (Sb/Sc). The 
primary, offset ~ 3" ESE from the SCUBA position, has a 
robust radio and 24/im detection, whereas the secondary, 



offset ~ 5" ESE has a robust 24/im detection but no radio 
detection. 

Lock850.064: There are three bright sources detected at 
IRAC wavelengths within 10" of the SCUBA position. The 
primary counterpart is the only one with significant 24/im 
flux although it is non-robust being offset by nearly 10" to 
the SSE. The source is detected in the radio (also 
non-robust) and is mildly blended with a neighbour of 
similar brightness at IRAC wavelengths. 
Lock850.073: The primary and secondary counterparts 
have consistent photometric redshifts and are well fit by 
late type SEDs. The primary shows an irregular optical 
morphology and has robust radio and 24/im detections. 
The secondary has a similar offset to the primary from the 
SCUBA position (~ 2.6"), has a robust radio detection but 
is not detected at 24/im. 

Lock850.077: This source has two possible counterparts. 
The primary gains its status by virtue of being close to the 
SCUBA position (~ 1.5" to the E) thereby having robust 
radio and 24/im detections. The brighter secondary (three 
times the flux in the radio and at 24/im) is offset 
approximately 6" to the ENE of the SCUBA position and 
the optical morphology clearly shows two very close (< 1") 
objects of equal flux and size. The 3" aperture we use 
completely encompasses both of these objects belonging to 
the secondary and therefore inaccurate deblending may 
explain the very poor SED fit. 

Lock850.083: The only plausible counterpart for this 
source is a nearby spiral galaxy with a photometric 
redshift of z = 0.22 offset ~ 8" W of the SCUBA position. 
The source has a robust 24/im detection but is not 
detected in the radio. This object is best fit with a type Sc 
SED, perfectly consistent with its optical morphology. 
Lock850.100: There are two very close potential 
counterparts for this source. The primary counterpart is 
the brighter of the two, with a photometric redshift of 
2 = 1.4 and having photometry consistent with a starburst 
galaxy. The fainter source is best fit with an Im SED at 
z — 0.4. The IRAC photometry is omitted for this object 
as the primary counterpart cannot be deblended with its 
neighbour due to the IRAC PSF. 



APPENDIX C: PHOTOMETRIC DATA 

Table IC II lists the optical to mid-IR photometry for all 
SHADES source counterparts discussed in this paper. All 
magnitudes are in the AB system and were extracted using 
an effective 3" diameter aperture (see Section[2]for more de- 
tails). All photometric errors were combined in quadrature 
with an error of 0.05 to account for uncertainties in zero 
point. 
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Table CI. The optical to mid-IR photometry for all SHADES source counterparts. Magnitudes are in the AB system and correspond 
to an effective 3" diameter aperture (see Section [2}. A numerical value of 99 corresponds to a non-detection, whereas a hyphen indicates 
that photometry could not be extracted. Errors include the zero point uncertainty in each waveband. 
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